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ABSTRACT 


Observations by the airborne X-band Doppler radar (EDOP) and the NCAR S-band 
polarimetric (S-Pol) radar from two field experiments are used to evaluate the surface reference 
technique (SRT) for measuring the path integrated attenuation (PIA) and to study attenuation in 
deep convective storms. The EDOP, flying at an altitude of 20 km, uses a nadir beam and a 
forward pointing beam. It is found that over land, the surface scattering cross-section is highly 
variable at nadir incidence but relatively stable at forward incidence. It is concluded that 
measurement by the forward beam provides a viable technique for measuring PIA using the 
SRT. Vertical profiles of peak attenuation coefficient are derived in two deep convective storms 
by the dual-w 7 avelength method. Using the measured Doppler velocity, the reflectivities at the 
two wavelengths, the differential reflectivity and the estimated attenuation coefficients, it is 
shown that: supercooled drops and (dry) ice particles probably co-existed above the melting level 
in regions of updraft, that water-coated partially melted ice particles probably contributed to high 
attenuation below the melting level. 
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1. Introduction 


Airborne and satellite-borne radars typically operate at wavelengths less than 3 cm to 
reduce the overall size and weight of the payload, and to obtain adequate spatial resolution. The 
most notable satellite-borne radar is the Tropical Rain Measuring Mission Precipitation Radar 
(TRMM-PR) which operates at a wavelength of 2.17 cm (Kummerow et al. 1998). The TRMM- 
PR is used to estimate rainfall, R, from the radar reflectivity, Z, by means of empirical Z-R 
equations. However, at these shorter wavelengths, the microwave radiation suffers significant 
attenuation in passing through precipitating storms. The attenuation is significant in rain and can 
be very pronounced in the presence of wet ice particles, such as, melting or water-coated hail or 
graupel. It is important to correct for the attenuation in order to estimate the rainfall more 
accurately. The study and observation of attenuation can help us to improve the algorithms for 
attenuation correction, and thereby improve the estimation of precipitation by the PR. For a 
given wavelength and polarization state, the microwave attenuation depends upon the size, 
concentration, shape, orientation and composition of the hydrometeors (Battan 1973). However, 
in situ microphysics are difficult to obtain in high reflectivity, highly attenuating core. Therefore, 
the measurement of attenuation can also help us to understand the microphysics of precipitation. 

Direct measurement of attenuation is difficult. Three methods have been used for both 
ground-based and airborne radars. The dual wavelength method measures reflectivity at the 
attenuating wavelength (e.g., 3 cm) and at a non-attenuating wavelength (e.g., 10 cm) 
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simultaneously in a common volume. The difference of the two reflectivity factors 1 gives the 
two-way integrated attenuation between the radar and the common volume along the path of the 
attenuating beam; this is usually referred to as the path -integrated attenuation (PIA). This method 
assumes that the wavelength dependence of the back-scattering cross-sections of the scatterers is 
accounted for; this can be done if the particles are in the Rayleigh scattering region (Eccles and 
Mueller 1971; Meneghini et al. 1994; Bolen and Chandrasekar 2000). The second method, called 
the dual-radar method, utilizes two-spaced radars operating at the attenuating wavelength to map 
the reflectivity field from two different aspects. The two reflectivity fields can be used to obtain 
the specific attenuation. This method has been used for airborne radar (Testud and Amayenc 
1989) and for ground-based radar (Srivastava and Tian 1996; Tian and Srivastava 1997). The 
dual-radar method does not require that the scatterers be in the Rayleigh range. The third method 
involves the use of a reference target of known reflectivity. When a surface target is used as the 
reference, the method is called the surface reference technique (SRT) (e.g., Meneghini et al. 
1983). In this method, the ‘reference’ radar cross-section of the surface is first determined in the 
absence of attenuation. In practice this means measuring the reflectivity of the ground in a 
precipitation-free area in close proximity to the precipitating area. When measurements of the 
ground cross-section are made through precipitation any decrease from the reference cross- 
section is attributed to two-way PIA between the radar and the surface. The SRT is used in 
processing TRMM-PR data (Iguchi et al. 2000). 

1 The equivalent reflectivity factor will hereafter be referred to simply as reflectivity for the sake of brevity. 
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The success and accuracy of the SRT depends upon the stability of the radar cross-section 
of the surface. A number of studies have been devoted to study the dependence of the surface 
reflectivity on the angle of incidence, especially over the oceans, which give a relatively constant 
microwave reflection. Limited studies of the SRT over land have been reported to date 
(Meneghini et. al. 2000). In general, this is because the radar cross-section of land surface can be 
highly variable particularly at near-nadir incidence angles. A second concern is the change in the 
surface cross-section between raining and non-raining conditions due to wetting and other 
surface changes accompanying precipitation. 

Several field experiments were conducted in support of TRMM. Of interest to us is the 
Texas and Florida Under Flights Experiment-B (TEFLUN-B) held near Melbourne, Florida in 
August and September 1998 and the TRMM Land-Biosphere Atmosphere (TRMM-LBA) held in 
southwestern Amazon in January and February 1999. 2 A number of universities and research 
laboratories participated in the experiments and deployed a variety of instruments. 

In this paper, we present selected instances of high PIA observed over land at high 
incidence angle using NASA’s ER-2 high altitude aircraft. The observations were obtained with 
a radar mounted on the ER-2 in Brazil and Florida dunng TRMM field campaigns. We compare 
the observations by the airborne X-band ER-2 Doppler Radar (EDOP) with nearly simultaneous 
and collocated observations obtained by a S-band ground-based radar (the National Center for 
Atmospheric Research, S-Pol radar) that is not subject to significant attenuation. The 



characteristics of the radars we use in this study are listed in Table 1. We shall compare the 
attenuation derived by the SRT from EDOP with attenuation derived by the dual wavelength 
method using EDOP and S-Pol described in Section 3 and 4. We shall interpret the measured 
attenuation in terms of microphysical properties of the precipitation in Section 5. 

2. Attenuation Observed by EDOP Using SRT 

a. SRT method 

We briefly review the method of obtaining PIA by the SRT and present an example of a 
storm that shows significant PIA. The observed return power from the surface, P s , is related to 


its normalized radar cross-section, G by: 
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for nadir incidence and 
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for forward incidence (Kozu, 1995), where C 5 and C sl are constants depending upon the radar 


characteristics and the angle of incidence. In the presence of precipitation, the apparent 
normalized radar cross-section of the surface, G R , is related to the normalized radar cross- 

section in the absence of precipitation, G° NR , by: 


2 http://www.eosdata.gsfc.nasa.gov/CAMPAIGN.DOC/TRMM-FE 
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where A is the two-way path -integrated attenuation and <7^ and crJ^j are in dB. Because O NR is 
not measurable, it is replaced by a reference value, an average of C7 NR along the flight 

track for an adjacent rain-free area. Using this reference value, the estimated PIA is: 



Because the reference value can be different from the true surface radar cross-section, A can 
differ from the true PIA, A. Differences between an< ^ ^ nr ma Y due viability 

of the surface, and changes of surface scattering due to changes in wind, vegetation, soil 
moisture, etc (Ulaby 1982). We need to know the magnitude of this variability in order to assess 
limitations in the measurement of the PIA. 
b. Application of SRT to EDOP data 

The Doppler radar (EDOP) mounted on the NASA’s ER-2 aircraft operates at 3-cm 
wavelength. It has two fixed antennas, one pointing at nadir and the second pointing 33.8° ahead 
of nadir. The antennas are identical with a beam width of 3 (circular) defining a footprint at 
surface of about 1 km at nadir assuming ER-2 flies 20 km above the surface. The ER-2 ground 
speed is about 210 m s" 1 and the integration time for the data processing is 0.5 s. These two 
values imply that the surface cross-section is over sampled with one sample being measured 
every 100 m along the flight track, and 10 samples are obtained over one beam w'idth. The range 
resolution of the radar is 37.5 m. Additional details about the radar and the data processing may 
be found in Heymsfield et al. (1996). 
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Figure la shows a histogram of the surface reflectivity observed by EDOP from the nadir 
and forward antennas for a cloud and precipitation free region over the ocean off the Gulf coast 
of Florida. The data are typical in that surface reflectivity of the ocean surface at nadir incidence 
has small fluctuations of + 1 dB while the echo at 33.8° incidence is more variable having a 
standard deviation of about 2 dB. The situation is very different over land (Fig. lb). At nadir 
incidence the reflectivity is highly variable having a standard deviation of about 4 dB. At 
forward incidence, however, the surface echo over land shows variability similar to that over the 
ocean with a standard deviation of about 2 dB. The importance of distributions such as those in 
Fig. 1 to the SRT is that the variability of the surface return limits the minimum PIA that can be 
measured and also gives an estimate of error in measured PIA due to surface variability. We see 
that over land, the SRT is subject to larger errors at nadir incidence, while at forward incidence 
PIA can be measured with the same accuracy as over the oceans. It should be noted that over 
land cr° can vary significantly depending upon terrain and vegetation (e.g., Ulaby, 1982), angle 
of incidence and radar characteristics such as side lobes. The terrain types observed during the 
observations presented here are mostly flat areas along the coast and the Amazon rain forest. We 
may conclude that in those regions the variability of < 7 ° over land at 33.8 incidence is about +2 
dB. Therefore, PIA greater than about 2 dB can be measured over land at 33.8° incidence angle. 
Caution should be used in applying this result to other types of terrain. 

It is emphasized that it is not necessary to measure the absolute value of the surface 
return since the SRT is a differential technique. It is important, however, that the magnitude of 
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the surface cross-section does not change as a consequence of wetting of the surface by the rain, 
since the SRT will attribute any such change to PIA. In the cases studied so far with the EDOP 
radar, no significant change in the surface cross-section has been observed at the transition 
between rainy and clear conditions for non-nadir incidence 

An example of the observation by EDOP is shown in Fig. 2 from an East- West flight 
over land near Ji Parana, Brazil on 12 February 1999. Figures 2a and 2b show the reflectivities 
measured by the nadir and forward antennas; Figure 2c shows the surface cross-section, <7 , 
measured by the nadir (dotted) and forward-pointing (solid) antennas. The flight track covers a 
convective cell embedded in stratiform rain. The reflectivity measured by the EDOP nadir 
antenna shows that, in the cell, the 10 dBZ contour reaches a height of 15 km and the peak 
reflectivity is 50 dBZ at a height of 4 km at a distance of about 42 km. 

A note about the terminology and coordinate system may be helpful here. In Fig. 2 and 
subsequently, the reflectivities measured by the EDOP forward and nadir-pointing beams will be 
denoted as Z (/ and Z xn respectively. Z s will represent the reflectivity measured by the S-Pol 

radar. All altitudes are above ground level (AGL) unless stated otherwise. The x-coordinate in 
the vertical sections is the distance along the flight path from some arbitrary origin; this is the 
same as the distance along the surface. It is to be noted that the forward reflectivity is measured 
when the radar is ahead of the x-coordinate shown in the figure; for the surface observation, it is 
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ahead by 20km x tan (33.8°), or 13.4 km. The PIA inferred from the forward beam is along a 

slant path that terminates at the surface at the x-coordinate shown. 

Returning to Fig. 2a, the nadir pointing antenna shows a low reflectivity region near the 
surface situated below a region of high reflectivity at a distance of about 40 km. The nadir 
surface cross-section, (7° (Fig. 2c - dotted line) shows that this low reflectivity is probably due to 
attenuation; however, it is difficult to interpret the decrease in the cross-section in terms of PIA 
with confidence because of the large fluctuations in the background surface cross-section at nadir 
incidence. The storm reflectivity measured by the forward-pointing antenna (Fig. 2b) is similar 
to the nadir reflectivity except that a lower reflectivity exists at about 45 km distance. The 
relative displacement of the minimum of the surface cross-sections for the forward and nadir 
beams is related to the beam orientations; the forward beam suffers greater attenuation along a 
slant path that intersects the ground ahead of the nadir beam. The integrated attenuation in the 
forward beam is greater probably because the slant paths are longer. In contrast to the nadir 
beam, the background <J° (regions of little or no precipitation) measured by the forward beam 
(Fig. 2c - solid line) is stable. The forward minimum <7° of -23 dB at 45 km, compared with the 
background (7° of -8 dB, yields a two-way PIA of 15 dB. It is encouraging to note that 
notwithstanding the large fluctuations in the nadir (7°, the structure of the dip in the nadir ( 7 ° is 
very similar to that of the forward (7°. The double dip in the nadir <7° seems to be correlated 
with the storm structure at high levels. 
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During the two TRMM field campaigns, the ER-2 flew over numerous deep convective 
storms. The cloud tops, defined by the 10 dBZ contour, reached up to about 13-18 km; many of 
these storms also showed large PIA. Table 2 lists eight such cases of storms over land. Most 
cases have rather intense convection as indicated by the fact that the 40 dBZ contour reached up 
tol5 km and the PIA exceeded about 20 dB; such PIA is large for the two-way path of about 20 
km. In the following, we present two cases from the table. One of the cases is from Florida and 
the other is from Brazil. These cases were selected because the S-Pol radar also observed these 
storms nearly simultaneously. S-Pol and EDOP can be used together to provide an independent 
estimate of attenuation using a dual-wavelength method. 

3. Dual- Wavelength Method and S-Pol Data Processing 
a. Dual-wavelength method and method of Using k-Z s relation 

The dual wavelength method takes the difference of the reflectivities observed by the S- 
Pol ( Z s , dB) and the EDOP ( Z x , dB) radars. In general, we can write: 

Z s -Z x =2l r 0 k{r')dr' + S, (4) 

where <5 is the difference in reflectivity due to the departure from Rayleigh scattering. For 
particles much smaller than the wavelength, <5 = 0. Figure 3 shows <5 for a mono-disperse 
distribution of spherical water drops of different diameters. We see that <5 can be taken as zero if 
the drop diameter is less than about 2.5 mm. A minimum 5 of about -3.5 dB occurs for drops of 
diameter about 6-8 mm. 
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Since 10 cm wavelength radiation is subject to little or no attenuation, an expected PIA at 3.2 
cm may be calculated from the power-law empirical k- Z s equation: 

it = aZ b s ,a = 2.9 x 10 _4 ,£> = 0.72 (5) 

by evaluating the integral of 2aZ b along the path of the EDOP beam. In equation (5), k is in dB 
km ' and Z s is in The values of a and b are taken from Battan (1973). Equation (o) is 

applicable for the modified Marshall-Palmer distribution, a wavelength of 3.2 cm and spherical 
raindrops having a temperature of 0 T. It may be noted that the PIAs or attenuations estimated 
from the above two methods, involving additional data from the S-Pol radar, are independent of 
the PIA deduced from the SRT. 
b. S-Pol data processing 

The NCAR S-Pol radar is transportable ground-based dual-polarization radar operated at 
10 cm wavelength. The beamwidth of the S-Pol antenna is 0.91 (circular) and its range 
resolution is 150 m. Besides Doppler parameters, S-Pol measures reflectivity and polarization 
quantities, such as, the differential reflectivity (ZDR) and the linear depolarization ratio (LDR). 
Further details for S-Pol may be found at the website. 3 

S-Pol radar observations were processed as follows. Data from the volume scan nearest in 
time and space to the ER-2 overflight was interpolated onto a grid coincident with the vertical 
plane mapped out by the ER-2 radar. The EDOP has generally higher resolution than S-Pol. For 

3 http://www.atd.ucar.edu. 
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example, at 50 km distance from the S-Pol, the approximate distance of S-Pol from the storm 
cases to be discussed, the resolution of the S-Pol beam is about 0.87 km in the cross-beam 
direction and 0.15 km in the along-beam direction; for the EDOP radar, the cross-beam 
resolution is about 0.76 km at 5 km height and the along-beam resolution is a constant 0.0375 
km. Spatial averaging of the data was performed in order to make their resolutions 
approximately comparable. Weighted arithmetic averaging of the reflectivity (mm 6 m 3 ) was 

performed and the dBZ value was calculated from that average. Spatial averaging of the data was 
performed using the Cressman filter, that is, a weighted average of the data was calculated using 
the weighting function ( R 2 - r 2 )/(R 2 + r 2 ) where R is a cutoff distance and r is the distance 

between the observation point and the point to which the interpolated value is assigned. Zero 
weight is assigned to observation points for which r>R. R was selected such that the 
resolutions of the interpolated values from the different observation platforms were nearly equal, 
it depended upon range from the radar and the radar characteristics. Such averaging should make 
the data from different platforms more comparable and reduce errors in estimating PIA or 
attenuation by the dual wavelength method. Still, we can expect errors in regions of high 
gradients and in the vicinity of storm boundaries. This is because of incomplete and non-uniform 
beam-filling, especially near storm boundaries. Errors can also be due to temporal displacements 
in data and such displacements are less than 2 mins; no attempt was made for time interpolations. 


12 



In the next section, we present data for the two storms and attempt to validate the 
attenuation inferred from the SRT. The attenuation is interpreted using a) dual-wavelength in 
terms of the microphysical structure of the storms and b) S-Pol multi parameter data. 

4. Case Studies of Large PIA 

a. A Storm on 15 August, 1998 in Florida 

On 15 August 1998, the ER-2 flew over a 60 km long N-S line of vigorous convective 
cells during 2223-2230 UTC. The merging of the East and West Coast sea breezes triggered the 
cells. Nearly simultaneous and co-located observations by the S-Pol were used to construct Fig. 
4, which shows (a) reflectivity, Z, (b) ZDR, and (c) LDR during 2224-2227 UTC at 8 km height. 
At this height, the maximum reflectivity is about 58 dBZ (indicated by +) about 1 km west of the 
ER-2 flight line whose projection is shown by the solid line. Coincident with the high reflectivity 
core, the ZDR is near zero (Fig. 4b) and the LDR is greater than -18 dB (Fig. 4c); this generally 

indicates the presence of wet hail (Doviak and Zmic 1993). 

Figure 5 shows vertical sections of the reflectivity, Z s , constructed from S-Pol radar data 

along the flight track (Fig. 5a), the reflectivity, Z xf , measured by the forward-pointing beam of 

the EDOP radar (Fig. 5b), the PIA deduced by the SRT (Fig. 5c, solid line), and the PIAs 
deduced by the two methods outlined above (Fig. 5c, dotted and dashed lines). We see that the 0- 
10 dBZ Z s contour reaches a height of about 15 km and the 40-50 dBZ contour reaches a height 
of about 14 km; these indicate an intense storm. The maximum reflectivities seen by S-Pol and 
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EDOP are 55 and 53 dBZ respectively (“+“in Figs.5a and 5b). Such high reflectivities suggest 
that hail may have been present in the storm. As pointed out above, wet hail was probably 
present west of the flight path. 

The maximum PIA measured by the SRT along the forward beam is 25 dB at a distance 
of about 72 km (Fig. 5c). The PIA obtained by integrating the empirical k — Z s equation (5) 
along the forward-pointing beam from storm top to 1 km height AGL is shown by the dotted 
line; its peak value is about 30 dB. [Reliable Z s were not available below 1 km height because of 
radar scan limitations and ground clutter.] The peak PIA from k — Z s relation is about 5 dB 
greater than that from the SRT. This discrepancy is greater than the estimated uncertainty (about 
2 dB) in the SRT. Moreover, the contribution of attenuation below 1 km height will make the 
discrepancy even worse. Therefore, we can probably attribute the difference between the PIAs 
deduced by the SRT and from the k-Z s relation to inadequacies in the empirical k~Z s 
equation, and presence of hydrometeors with characteristics other than those on which the 
empirical k-Z s equation is based. 

The PIA deduced from the difference, Z s - Z xj , along the forward beam from flight 

altitude down to a height of 1 km (dashed line) generally parallels the other two PIA curves. Its 
peak value, 28 dB, is in good agreement with the peak PIA deduced from the SRT. This 
agreement is quite remarkable considering that the data came from two radars in different 
locations and having different resolution volumes. A noticeable difference between the dual- 
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wavelength PIA and the other PIAs is the rather pronounced dip in the dual-wavelength PIA at a 
distance of about 74 km. This difference may be due to the interpolations performed on the data 
that are more prone to produce biased results near storm boundaries as discussed earlier. Another 
possibility is that the dip is due to dual-wavelength effects, that is, non-Rayleigh scattering by 
large hail or raindrops (for the latter, see Fig. 3). This is to some extent supported by the fact 
that, along the slant path of the radar beam associated with the point where the dip occurs, there 
is a pocket of high downward Doppler velocities of magnitude about 15-20 ms as will be seen 
later (Fig. 15b, pocket at about 71 km at height of 3-4 km) suggesting large particles and / or 
strong downdrafts. We favor the notion of large raindrops, rather than hail, because the ZDR 
values are positive in this region (about 1.5-2. 5 dB) and show an upward bump (Fig. 15a also 
to be discussed later). However, this suggestion should be regarded as tentative in view of the 
difficulties associated with processing two such diverse sets of radar data. 

It is also possible to estimate range-resolved PIA by taking the differences of Z 5 and Z ^ 

at different heights. This procedure would yield the attenuation along slant paths. From a 
microphysical standpoint, vertically resolved PIAs are of greater interest, especially in deep 
convective clouds. Therefore, we shall concentrate on range-resolved attenuation inferred from 
Z s and Z xn , the reflectivity observed by the nadir beam of the EDOP. 

Before presenting the range-resolved attenuation, we compared the reflectivities 
measured by the S-Pol and EDOP radars because the accuracy of the dual-wavelength method of 
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PIA estimation depends upon the accuracy of the radar calibration besides the assumption of 
Rayleigh scattering. Figure 6 is a scatter plot of Z s vs Z xn for the storm in Fig. 5. Except for a 
relatively small fraction of the points, Z, and Z xn scatter around the one-to-one line. The 
standard errors of Z t and Z xn are about 1-2 dB, implying a standard error in Z s —Z xn of about 
1.5- 3 dB. This is confirmed by the scatter plots. Therefore, we should be able to measure PIAs 
greater than 2-3 dB by the dual wavelength method with a 2-3 dB accuracy. It should be noted 
that the accuracy of the range-resolved attenuation does not depend upon the absolute accuracy 
in the measurement of the reflectivities. 

Figure 7a is a plot of the S-Pol reflectivities, Z s , along the paths of the nadir beam at the 
heights of 1, 2, 3, 4, 5, and 6 km. Figure 7b shows dual wavelength PIAs for the nadir beam from 
storm top to those heights. Figure 7c shows the PIA from integration of the empirical k — Z s 
equation along the nadir paths. In Fig. 7a ,we have plotted only reflectivities greater than 45 dBZ 
because lower reflectivities contribute insignificantly to the PIA; this can be seen by comparing 
Fig. 7a with Figs. 7b and 7c. The PIA curves in Figs. 7b and 7c have similar shapes but differ in 
magnitude. The two sets of curves can be made to agree better by ‘tuning’ constants a and b in 
the (5) which depend on the drop size distribution. The main peaks of the PIA in Fig. 7b occur at 
a distance of about 69 km in the region of highest reflectivity. A second peak occurs at about 75 
km in the region of sharp gradient of reflectivity. The PIA in the latter region may not be reliable 
because of the reasons mentioned earlier. It is encouraging to note that the curves for the higher 
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heights are generally below the curves for the lower heights. While this must necessarily be so in 
the case of Fig. 7c, it is not necessary in case of Fig. 7b. In the case of Fig. 7b, this tendency can 
be reversed in the presence of non-Rayleigb scatterers, such as large hail, or significant signal 
fluctuations. The fact that the curves for lower heights do lie above the curves for higher heights 
in Fig. 7b strongly suggests that large hail was not present. It may be mentioned that this does 
not contradict our earlier discussion (Fig. 4) because there we found evidence for hail about 1 km 
to the left of the plane of observation by the EDOP radar. 
b. A Storm on 10 February, 1999 in Brazil 

The EDOP flew over a line of convective storms in Brazil during 1811-1816 UTC,10 
February 1999. The southern most cell weilTlocated about 50 km north of the S-Pol radar. Figure 
8 shows S-Pol reflectivity at 3 km height; it was constructed from the 1810-1815 UTC S-Pol 
volume scan. 

Figure 9a, similar to Fig. 5a, indicates that 0-10 dBZ and the 40-50 dBZ contours reach 
heights of about 15 km and 7 km, respectively. The maximum reflectivity (“+” in Fig. 9a) is 52 
dBZ, lower than the Florida storm discussed previously. The maximum reflectivity from forward 
beam (Fig. 9b) is only 50 dBZ. We also note that the shape of the lower reflectivity region (Fig. 
9b, 30-32 km) is aligned in the direction of the forward beam, a signature of attenuation (Atlas 
and Banks, 1951). The two-way PIA in the direction of the forward-beam derived by the SRT 
(Fig. 9c, solid line) shows a maximum PIA of 29 dB at a distance of about 32 km. The PIA along 


17 



the forward path from storm top to 1 km height from the empirical k Z s equation (Fig. 9c, 
dashed line) is similar in shape to the SRT curve, but is substantially lower; a peak value of 
about 12 dB at 32 km distance versus 29 dB for SRT. The PIA to 1 km height along the forward 
beam deduced by the dual wavelength method (Fig. 9c, dotted line) also parallels the SRT but its 
peak value of about 19 dB is much smaller than the PIA from the SRT. This peak would have 
been closer to the maximum SRT value of 29 dB if the attenuation below 1 km height could have 
been included, but would still fall short of the maximum SRT value. Other discrepancies 
between the dual wavelength PIA and the SRT PIA occur mainly in regions of high gradients 

and at storm edges for reasons already noted. 

Figure 10 is a scatter plot of Z s against Z xn , similar to Fig. 6. The X-band and S-band 

reflectivities have more scatter about the one-to-one lines than that the Florida storm (Fig. 6). 
One probable reason is that for the previous Florida case, an RHI scan was available in close 
proximity to the EDOP flight plane and time. The Z s shown was essentially obtained by 
interpolating the Z s values from the RHI grid to the grid of the EDOP observations. In the case 
of the Brazil storm, however, a complete volume scan with PPIs at different elevation angles had 
to be used to construct the vertical plane. As stated before, the volume scan took about 5 min. In 
the EDOP observations, there is a systematic progression of time with horizontal distance during 
data acquisition leading to a systematic difference between the times at which the S-Pol and 
EDOP observations are acquired. No attempt was made to correct for differences in time. We 
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believe that the greater time difference between the two radar observations and the a lower 
density of observation points contributed to the greater scatter in the Brazil case. 

Figure 11 for the Brazil storm is similar to Fig. 7 for the Florida storm. We again note 
that the PIA curves derived from the empirical k-Z s equation (Fig. 1 lc) are similar to the dual- 
wavelength PIA curves (Fig. lib) but they are considerably smaller in magnitude. As mentioned 
previously, better agreement could be obtained by tuning the coefficient and exponent in the 
empirical k-Z s . 

The peak two-way PIAs and range-resolved attenuation along a vertical path (vertical 
solid line in Fig. 7b and Fig. 1 lb) for Florida and Brazil storm are summarized in Table 3 and 4 
for later discussion. The Z and ZDR in the tables were obtained from Fig. 16 and 18 to be 
discussed later in section 5b. In Table 4, there is a negative attenuation in the height interval, 1.5- 
2.0 km which is larger than the 1.5 dB uncertainty estimated in the dual-wavelength reflectivity 
ratio. In deriving the range-resolved attenuation, we take the difference of the reflectivity ratio at 
the two heights. Also, as noted earlier, the observations for Brazil storm are probably subject to 
greater processing errors than the Florida storm. Thus the negative attenuation may be apparent, 
being the results of various sources of error. Another possibility is that non-Rayleigh scatterers 
were present in significant concentration in this height interval. 

5. Inferred Microphysics from Theory and Observations 
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In both the Florida and Brazil storms, the attenuation between the storm top and 5 km 
height, estimated by the dual-wavelength method, was 1 dB or less (Table 3 and 4). The storm 
top ranged from 10-15 km AGL, ensures that ice particles were present in the storm. The small 
attenuation above 5 km means that this region was dominated by ice particles and that wet ice 
particles which attenuate strongly were not present in significant concentrations. As mentioned 
before, LDR from S-Pol showed evidence of wet hail at an altitude of 8 km in the Florida storm, 
but this was to the West of the flight path of EDOP. In that region, the polarimetric data (not 
shown) also showed evidence of wet ice particles below the melting level. Furthermore, large 
hail greater than about 1 cm in diameter was also probably not present in concentrations large 
enough to cause significant attenuation because it would have been manifested by its Mie 
scattering effects at the shorter wavelength. Interpretation of the attenuation coefficients inferred 
for the lower levels is more problematical because raindrops and mixed phase particles were 
likely to have been present below the melting level. 

In this section, we attempt to interpret the peak attenuation estimated above in terms of 
the hydrometeors responsible for the attenuation. For this purpose, we shall primarily use the 
attenuations obtained by the dual-wavelength method, since this method yields range-resolved 
attenuation coefficients. Because of the great complexity and diversity of the possible 
distributions of hydrometeors in deep convective storms, conclusions from our limited data set 
must be considered to be qualitative and preliminary, since there is no direct validation. 


20 



In order to better understand the process below the freezing level, we also use ZDR 
observed by the S-Pol radar and the Doppler velocity, V n , observed by the nadir beam of the 
EDOP radar. Further, we shall make use of theoretical relationships between reflectivity, 
attenuation coefficient and ZDR for Gamma function distribution of raindrop sizes. These 
relationships provide a ‘base reference’ case for interpreting the data. We begin with the 
theoretical relationships. 
a. Theoretical Relationship 


1) Attenuation, Reflectivity and Median Volume Diameter: 

We assume a gamma rain drop size distribution (Ulbrich, 1983): 

N(D) = N 0 D t* exp(-AD) (6) 

where N(D)AD is the concentration of drops of diameter D to D + A.D and N 0 , jA and A are 


parameters of the distribution. The exponential distribution is a special case of the (6) (p 0). 
To facilitate analytical manipulations, we assume that the diameter ranges from zero to infinity 
in (6); qualitatively similar results are obtained when a non-zero minimum diameter and a finite 
maximum diameter are considered. The median volume diameter, D 0 , of the above distribution 

is given approximately by: 


3.67 + p 


(7). 


The reflectivity factor is given by: 
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( 8 ), 


Z = c, 


A 


lKZ±ifl 


A 


The extinction cross-section of a spherical raindrop can be approximated by a power-law 


equation (Atlas and Ulbrich, 1974): 


Q=aD n 


(9) 


where the coefficient a and the exponent n depend upon the wavelength and temperature. Using 
the above equation, the attenuation coefficient is given by 


k = c, 


a T N, 


A" JLA m 


r(« + n + 1 ) 


A 


( 10 ) 


From (7), (8) and (10) we can deduce: 

z c z r(7 +■ jij i 


d: 


6 -n 


k c k + pi + l) ci (3.67 + fij 


6-n 


(ID 


In the above equations, c z and c k are numerical constants and the equations have been written 
such that terms in square brackets have rational dimensions (no fractional exponents of the 
length are involved). The values of the coefficients and the dimensions of the various terms are: 


[Z] = 


mm 


m 


,[&] = dBkm 1 


A 


= [cm], 


a 

A 7 


= [ cm2 l 77 = [ cm ‘ 4 l 


A M 


( 12 ) 


[Z7 0 ] = [cm],[lV 0 ^] = [cm- 4 ] 
[Q t ] = [cm 2 ],[D] = [cm] 



c. = 10 12 ,q =4.343 xlO 5 

The coefficient a and the exponent n in (9) are taken from Atlas and Ulbrich (1974, 

0 

Table 1). For 3 cm wavelength, a = 4.18 to 13.24 and n = 4.33 to 5.16 for 0 to 40 C 

From (1 1), we see that for given D 0 and jJ . , Z and k are proportional to each other with a 

proportionality factor that depends upon the temperature through a and n . Using (7), (8), and 
(10), we have prepared Fig. 12, which is a plot of k vs Z for four temperatures (0, 10, 18 and 

40 °C), three median volume diameter, Dq = 1, 2 and 3 mm, and three values of /r, namely, -2, 0 
and 4. Observed DSDs usually have \l in the range -2 to 3 and N 0 and fi are correlated. In 
preparing Fig. 12, it was not necessary to assume a value of N 0 because for a given u and Dq, Z 
determines N 0 (see (7) and (8) ) and k can then be calculated for any given temperature using 
(10). We see from (11), other factors being equal, increasing temperature generally gives greater 
k because a and n in (9) increase with temperature; however, this trend is reversed when small 
drops predominate, that is, for small Dq. We see that, for a given Z, the median volume 
diameter, Dq, has a pronounced effect on k, a decrease of Dq results in a significant increase in 
k. The inferred attenuation coefficients, listed in Tables 3 and 4, are also plotted on Fig. 12 and 
will be discussed later. In this and subsequent figures, no attempt has been made to depict error 
bars for the probable uncertainties in the inferred k; the lengths of the bars represent the range of 
reflectivities in the height interval, the ranges having been obtained from Figs. 16 and 18 to be 
discussed later. 
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2) Median Volume Diameter and Differentia] Reflectivity 

The differential reflectivity (ZDR) is the ratio of the reflectivities at horizontal and 
vertical polarizations. Raindrop shapes are well approximated by ellipsoids of revolution with 
their axis ratios being a function of the drop volume. Using equations given in Seliga and Bringi 
(1976), and the raindrop axis ratio given by Andsager et al. (1999, Eq. 1), we have calculated and 
plotted ZDR in Fig. 13 as a function of D 0 for the gamma DSD for four values of , u. Use of 
Andsager et al.'s formulation of drop shape is not crucial. Different formulation produce only 
small difference in the calculated ZDR and we are interested mainly in patterns of variation of 
ZDR and in diagnosing the presence of liquid drops. A maximum drop diameter of 6 mm was 
assumed in the calculations. The range of median volume diameters plotted is limited by the 
range of slopes, A , considered to be acceptable, namely, 16 to 40 cm' 1 . We may note that, unlike 
Z and k, both ZDR and D 0 are independent of N 0 , the absolute magnitude of the DSD, and 
depend only on its shape. On this figure, we have also indicated the ranges of the observed ZDRs 
for the different height intervals for the two storms, The ZDRs were taken from Figs 16 and 18 
presented later. Along the horizontal axis, the observed ZDRs have been positioned to intersect 
the curves for one or more values of . These observational points will also be discussed later. 

3) Reflectivity factor, Attenuation coefficient and Differential Reflectivity 

Both variables, Z and k, in Fig. 12 are proportional to N 0 [Eq. 8 and 9], while ZDR in 

Fig. 13 is independent of the magnitude of the DSD. Figure 12 can be used to estimate Z) 0 and 
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^ of gamma DSDs that are consistent with the observed Z and k. Figure 13 can be used to infer 
values of D 0 and fi of gamma DSDs that are consistent with the observed ZDRs. In the latter 
case, the magnitude of the distribution can then be inferred from the observed Z . However, Fig. 
13 makes use only of the ZDR from the S-Pol radar while Fig. 12 makes use of Z and k; the 
results from the two need not agree. To use all the three measurables, namely, Z, k and ZDR 
conveniently and consistently, we have constructed Fig. 14 which is a plot of Z/k vs ZDR for 

H - -2,0,2 and for temperature = 0, 10, 18, 40 C. Both Z/k and ZDR are independent of Nq. 
However, together they make use of the three measurables, namely Z and ZDR observed by the 
S-Pol radar at the 10 cm wavelength, and k inferred from the S-Pol and EDOP reflectivities. It 
may be noted that Fig. 14 is not independent of Figs. 12 and 13 but a convenient combination of 
the two figures that will facilitate our discussions. On Fig. 14, we have also plotted the 
observations for the Florida and Brazil storms; the Z/k and ZDR values are given in Tables 3 and 
4. 

b. Vertical Doppler velocity and ZDR Structure of the Storms 

Vertical sections of the reflectivities of the two storms have been presented in Fig. 5 and 
Fig. 9. We now present additional observations of the storms to facilitate the later discussion of 
the inferred attenuation coefficients in the context of the vertical structure of the storms. 

Figure 15 shows gray shade contour plots of (a) ZDR and (b) V n , the Doppler velocity 

measured by the nadir-pointing beam of the EDOP radar, for the Florida storm. Contours of Z s 



have been superposed, and the dashed line shows the locus of maximum V n . In Fig. 16, we show 
three vertical profiles of Z s , Z xn , V n and ZDR. The central plot is through the region of peak 

attenuation and the other two plots are 1 km on either side of it. 

Figure 15b shows that this storm has a sloping updraft. A maximum upward V n greater 

than 10 m s ' occurs at a height of about 11 km at 66 km distance. This implies vertical air 
velocity greater than 10ms’‘ in that region since V n is the resultant of the vertical air velocity 

and the reflectivity-weighted fall velocity of the scatterers. In the vertical profiles we note a 
regular progression of the region aloft of V n > 0 in going from the left panel to the right panel 
which is a reflection of the sloping updraft. We note that below the melting level the Z s is 
practically constant (central panel, Fig. 16) except in the lowest levels. We also note that below 
the melting level, the ZDR contours slope in the same sense as the updraft. The ZDR increases 
downwards from about 1.7 dB at the melting level to about 3 dB at the 1 km level (see central 
panel, Fig. 16; we are concentrating on the central panel because it is the region through which 
attenuation was derived earlier); an exception occurs in the 1-2 km height interval in which the 
ZDR shows a small decrease downwards. 

Figures 17 and 18 for the Brazil storm are respectively similar to Figs. 15 and 16 for the 
Florida storm. This storm has a structure similar to that of the Florida storm but it is weaker as 
indicated by the reflectivities. The V n contours again show a well-developed sloping updraft. 

Upward Doppler velocity greater than 10ms' occurs at about 10 km height and a distance of 33 
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km, implying an updraft greater than 10 ms in that region. We see again that below the melting 
level, the ZDR contours slope in the same sense as the V n contours aloft and Z s is again 
practically constant. The ZDRs below the melting level are smaller in this storm compared to the 
Florida storm implying that the drops in this storm were much smaller than the drops in the 
corresponding region of the Florida storm in keeping w'ith the weaker character of the Brazil 
storm. In contrast to the Florida storm, the vertical profile of ZDR (central panel) show's that the 
ZDR was practically constant with height below the melting level (excluding the lowest level 
where the data are questionable). 
c. Inferences from the observed Z, ZDR and k 
1) Florida Storm: 

From Fig. 12, we see that the observed k and Z for the 5-6 km (F5-6) height interval 
could be explained in terms of a gamma rain DSD with /x = -2 and a high D Q of about 3 mm 

(Fig. 12a). On the other hand, the observed ZDR in Fig. 13 is consistent not only with a rain DSD 
with = -2 and Z) 0 of about 1 mm, but also with DSDs with fx = 0 and D 0 in the range 1 .2 to 
1.6 mm, /X = 2 and D 0 in the range 1.4 to 1.9 mm, and // = 4 and D 0 in the range of about 1.9 
to 2.1 mm. None of the these DSD parameters are consistent with D 0 > 3 mm suggested by Fig. 
12. This is evident in Fig. 14 where here the observed Z/k and ZDR for F5-6 are not consistent 
with any gamma DSD considered. The main reason for this inability to explain the observations 
in terms of a gamma rain DSD is that the inferred k is too small for the observed Z, giving a 
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large value of Z/k. In view of the fact that this height interval is mostly above the melting level 
(4.8 km), we conclude that this region was dominated by dry ice particles causing little 
attenuation. However, the observed ZDRs above the melting level (0.3 to 1.7 dB in the 5-7 km 
height interval) imply that super-cooled drops must have been present in that region. 
Microphysically, this is supported by the presence of a strong updraft in this region. Thus, we 
may conclude that super-cooled drops and ice co-existed in this region with the liquid drops 
providing the small attenuation that w'as observed. 

The next lower height interval, F4-5 has its top near the melting level; therefore liquid 
drops and partially melted ice particles may be expected in this region. We see from Fig. 12 that 
the observed Z and k are consistent with a gamma DSD with fi= -2 and D 0 slightly greater than 
1 mm (Fig. 12a), n = 0 and D 0 slightly less than 2 mm (Fig. 12b), and a DSD with fX = 4 and a 
Dq somewhat greater than 2 mm (Fig. 12c). From Fig. 13, we see that the observed ZDR is 
consistent with DSDs with fl= -2 and D 0 -1 mm; fi = 0 and D 0 ~1.6 to 2.0 mm; n = 2 and 
D 0 ~ 2 to 2.4 mm; and ju = 4 and D 0 ~ 2.2 to 2.6 mm. However, Fig. 14 shows that the observed 
Z/k and ZDR are not consistent with any of the gamma DSDs considered. This is because the k 
is too high for a gamma distribution with the observed Z, giving a point that lies below the 
curves in Fig. 14. We conclude that this layer contained a mixture of water drops and water- 
coated ice particles formed by the partial melting of ice particles falling across the 0°C level. 
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Similar to F4-5, the observed points for F3-4 and Fl-2 lie below all the curves in Fig. 14. 
The attenuation is too large again to be due to a gamma DSD. Therefore, we conclude that these 
regions also contained partially melted ice particles that contributed to the high attenuation. It is 
important to note that if only some of the observations such as Z and ZDR had been considered, 
the observations would be falsely interpreted in terms of a gamma DSD. The additional 
information provided by the X-band observations, through the inferred attenuation, strongly 
suggests that mixed phase particles, rather than just water drops, were present in this region. 

Our hypothesis of melting particles is also supported by the systematic increase in ZDR 
with decreasing altitude as seen in the vertical section and profile of ZDR (Figs. 15 and 16) as 
well as Fig. 14. The observed increase in ZDR implies a substantial increase in the mean drop 
size with decreasing altitude. If only raindrops had been present, the observed increase in mean 
drop size could not be explained because the evolution of rain DSD by coalescence, which tends 
to increase the drop size, is very small and is counteracted by drop breakup. We believe that the 
observed increase in ZDR could be due to the formation of larger drops through the melting of 
larger ice particles with distance down. Flowever, this hypothesis needs to be tested through 
microphysical calculations of the melting of ice particles; this task is left for the future. 

We assumed that partially water-coated ice particles provided the observed high 
attenuation. This is supported by earlier work (Battan and Herman 1962). We have calculated the 
extinction cross-section of melting spherical ice particles in which the melt is assumed to form a 
concentric w'ater coat. This is obviously a simplification, because the coat may be eccentric, the 
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ice core may be not spherical, and the melted water may soak inside the particle, especially in the 
case of a particle of low bulk density. However, the simplified calculations should provide a 
qualitative guide. Fig. 19 shows the normalized extinction cross-section as a function of the 
fraction of mass melted for several values of the melted diameter and three values of bulk ice 
density. The extinction cross-section is normalized by the geometric cross-section of the melted 
particle. We see that the extinction cross-section of the partially melted particle can easily exceed 
that of the melted particle for certain sizes and melt fractions by a factor of several. The effect 
could be larger for deformed particles (Atlas et al. 1953). 

2) Brazil Storm 

The results for this storm are very similar to those for the Florida storm. From Fig. 12, we 
see that for the 4-5 km height interval (B4-5) which straddles the melting level at 4.2 km AGL, 
the observed Z and k are explainable by a gamma rain DSD with fJ. = -2 and a large D 0 of about 
3 mm (Fig. 12a). This high D 0 is inconsistent with the observed ZDR; according to Fig. 13, the 
ZDR could be due to gamma DSDs with (J. = -2 and D 0 — 1 mm, fi = 0 and Dq = 1.5 mm, ,u — 
2 and D 0 = 1.8 mm, and ^ = 4 and D 0 = 2 mm. This inconsistency between Figs. 12 and 13 is 
also evident in Fig. 14, which shows that the Z/k lies above all the points considered. We again 
attribute this to a small attenuation because of the predominance of dry ice particles above the 
melting level. However, the observed ZDR (~1.5 dB) again suggests that some super-cooled 
drops were present; these drops account for the small attenuation that was inferred. Again, the 
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presence of a strong updraft in this region supports the hypothesis of liquid drops above the 
freezing level. We note the interesting and perhaps significant fact that for both storms, the 
points representing the layers immediately above the respective melting levels occupy very 
similar positions on the diagram of Fig. 14. 

The observations for the lower layers, B3-4 and B2-3, are very' similar to the observations 
in the Florida storm below the melting level. Similar to the discussion to the Florida storm, for 
Figs. 12 and 13, two of the three parameters, Z, k and ZDR can be used to select the parameters 
of a gamma DSD that fits the observations. However, Fig. 14 again shows that the Z/k points he 
below all the gamma DSDs considered. This implies a k too high for the observed Z to be 
explained by a gamma DSD. As in the case of the Florida storm, we conclude that liquid drops 
and partially melted, water-coated, ice particles were responsible for the high attenuation. 

The ZDR was nearly constant below the melting level in contrast to the Florida storm, 
which showed a systematic and considerable increase with decreasing altitude. The ZDRs for the 
Brazil storm are also considerably smaller than for the Florida storm. This implies that the Brazil 
storm had smaller particles. This is consistent with the lower intensity of the Brazil storm. 

6. Summary and Concluding Remarks 

Observations by the EDOP (NASA’s ER-2 Doppler) radar, which operates at 3.2 cm 
wavelength and makes observations in a nadir and a 33 forward-pointing beam, show that over 
land the surface cross-section is highly variable at nadir incidence but is stable to within 1-2 dB 
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at the forward incidence. There was no evidence of changes in the surface cross section before 
entering and after exiting regions of heavy precipitation. Thus, measurement by the forward 
beam provides a viable method for measuring PIA (path integrated attenuation) by the SRT 
(surface reference technique) over land. Observations over land in a number of deep convective 
storms in Florida and Brazil yielded PIAs exceeding 20 dB between storm top and surface. We 
selected two storms, one in Brazil and one in Florida, for detailed study using the EDOP and S 
Pol (NCAR S-band Polarization) radar data. 

The PIA between storm top and ‘surface’ w-as also estimated by the dual-wavelength 
method using the S-Pol (S-band) and EDOP (X-band) reflectivities. In addition, the dual 
wavelength method was used to estimate range-resolved specific attenuation using the S-band 
radar data and the EDOP radar data from the nadir beam. Very good agreement was found 
between the forward beam PIA from the SRT and the dual wavelength method for the Florida 
storm; the agreement was not as good for the Brazil storm. The PIA was also calculated by 
integrating an empirical equation relating attenuation coefficient and reflectivity (S band) for rain 
DSDs (drop size distributions). It was found that the PIA so deduced and the PIA found by the 
SRT had very similar shapes though they differed in magnitude; the two PIAs could probably be 
made to agree quite well by ‘tuning’ the coefficient and exponent in the empirical attenuation- 
reflectivity equation. However, this tuning procedure has no physical basis as shown by the 
inferences from the dual-wavelength data and the observed differential reflectivities. 
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Both storms showed small attenuation above the melting level and significant ZDRs in 
layers a few kilometers thick lying just above the melting level. We concluded that this region 
was dominated by dry ice particles, but had coexisting supercooled drops in order to account for 
the observed ZDRs. The existence of supercooled drops is also supported by the Doppler 
velocity observed by the nadir beam of EDOP. The Doppler velocities showed strong sloping 
updrafts above the melting level. 

Below the melting level, rather large specific attenuations were deduced by the dual- 
wavelength method. Theoretical relationships between the attenuation coefficient, reflectivity 
and ZDR for a gamma rain DSD were formulated and discussed at length in an attempt to 
interpret the observations in terms of gamma rain DSD. It was possible to explain the 
observations in certain height intervals by means of a gamma rain DSD and deduce its 
parameters if only two of the three parameters were considered. However, such interpretation 
was found inconsistent with all three parameters mainly because the inferred attenuation, below 
the melting level, was too large for the observed Z. We also found a marked and systematic 
increase in ZDR with decreasing altitude, in the Florida storm, that can not be accounted for by 
the evolution of a rain DSD. It was suggested that the large attenuation was due to partially 
melted water-coated ice particles that present a larger extinction cross-section than the 
completely melted particle. Calculations for spherical water coated particle showed that the gain 
in extinction could be quite pronounced for ice particles of low bulk density. It is also known that 
deformed particles can cause even more extinction. The Florida storm was more intense than the 
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Brazil storm as suggested by its higher reflectivity. It contained bigger particles, as suggested by 
larger ZDRs; the larger ice particles required longer fall paths for complete melting and were 
probably responsible for the observed downward increase of ZDR. 

This study has shown the benefit of combining X-band airborne radar observation and S- 
band polarmetric measurements, i.e., a dual-wavelength method, for the study of microphysics of 
convective storms. Earlier attempts at using dual-wavelength method for ground-based 
measurements have used collocated radars, or common antennas with the matched beamwidths 
to overcome problems anticipated from lack of coincidence of radar resolution volumes. 
However, those methods usually yield the attenuation along horizontal paths. For deep 
convective storms, the structure along a vertical path is of paramount interest. Our study has 
shown that the combination of a down looking short wavelength (Doppler) radar and a ground- 
based S-band polarimetric radar can be used to infer plausible attenuation coefficients in deep 
convective storms. Valuable inferences about storm microphysics can be drawn by combining 
resulting dual-wavelength and ground-based polarimetric data. Such studies are important in 
order to improve the SRT approach over land for TRMM future satellite and to provide a firmer 
physical basis for its use. With proper design of scan strategies to obtain optimum space-time 
coincidence between ground-based and airborne radar data sets, it should be possible to obtain 
more accurate results. The interpretation of dual wavelength reflectivity, polarimetric, and 
Doppler velocity observations needs to be supported by further studies of the characteristics of 
melting particles, in particular, their radar and extinction cross-sections. 
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FIGURE LEGENDS 


Figure 1. Histogram of surface scattering cross-section over (a) ocean and (b) land at 
forward (solid) and nadir (dashed) incidence observed by the EDOP radar. 20 dB are 
added to forward surface cross-section for display purpose. 

Figure 2. Vertical cross-section of reflectivity at (a) nadir and (b) forward incidence, and 
(c) surface cross-section for forward (solid) and nadir (dashed) beams observed by EDOP 
radar on 12 February, 1999 in Ji Pama Brazil. 

Figure 3. Reflectivity at S-band (solid), X-band (short dashes) and their difference (long 
dashes) for mono-disperse distributions of spherical raindrops as a function of drop 
diameter. 


Figure 4. Constant altitude, 8 km agl, contour maps of (a) reflectivity, (b) differential 
reflectivity, ZDR, and (c) linear depolarization ratio, LDR in a storm in Florida on 15 
August 1998, 2224-2227 UTC reconstructed from S-Pol radar data. The location of 
maximum reflectivity is marked by + and the projection of the flight path of the airborne 


EDOP radar is shown by the line. 



Figure 5. Vertical cross-sections corresponding to Fig. 4 for the Florida storm, (a) 
reflectivity, Z s , from S-Pol radar, (b) forward beam reflectivity, Z^., from EDOP radar, 
and (c) path integrated attenuation, PIA. deduced using the surface reference technique 
(SRT, solid), integration of empirical equation 4 (dotted) and the dual wavelength method 
(dashed). .The 0°C level is indicated by dotted line. 

Figure 6. Scatter plots of S-band reflectivity, Z s , against X-band reflectivities at the nadir 
incidence for the Florida storm. The one-to-one lines are also shown. 

Figure 7. Plots for the Florida storm, for the heights (AGL) of 1.0, 2.0, 3.0, 4.0, 5.0 and 6 
km: (a) S-band reflectivity, Z 5 (b) difference of S-band and nadir beam X-band 
reflectivity, Z s - Z xn and (c) PIA from integration of empirical equation 4 along the nadir 
beam. Vertical line in the plot indicates the distance at which the maximum attenuation 
occurs. For details see text. 

Figure 8. Constant altitude, 3 km agl, contour map of reflectivity for a storm in Brazil on 
10 February 1999, 1810-1815 UTC, reconstructed from S-Pol radar data. The location of 



maximum reflectivity is marked by + and the projection of the flight path of the airborne 
EDOP radar is shown by the line. 

Fig. 9. Similar to Fig. 5 but for the Brazil storm. 

Fig. 10. Similar to Fig. 6 but for the Brazil storm. 

Fig. 11. Similar to Fig. 7 but for the Brazil storm. 

Fig. 12. Specific attenuation at 3 cm wavelength versus reflectivity at 10 cm wavelength 
for gamma function drop size distribution with shape parameters of (a) fJ. — — 2, (b) 
H = o, and (c) jX~ 4, for median volume diameters D 0 = 1, 2, 3 mm and temperatures of 

0, 10, 18 and 40°C. The inferred specific attenuations and reflectivity along the vertical 
line in Fig. 7 and 11 for the Florida and Brazil respectively are also plotted against 
reflectivity observed by the S-Pol radar with symbols Fx and Bx where x denote the layer 
interval in km. 


Fig. 13. Differential reflectivity, ZDR, against median volume diameter for fi - -2, 0,2, 4 
for the gamma drop size distribution. The range of median volume diameters plotted is 



limited by the range of slopes considered, namely, 16 to 40 cm' 1 . Fx and Bx specify the 
Florida and Brazil storms, respectively, where x denote the layer interval in km. For 
details, see text. 

Fig. 14. Plot of ZJk versus ZDR for gamma rain DSD for (a) n = -2, (b) [i = 0 and (c) 
[i = 2, for four temperatures. The observed ‘points’ are also depicted on the figures. Fx 
and Bx specify the Florida and Brazil storms, respectively, where x. denotes the layer 
interval in km. For details, see text. 

Fig. 15. Vertical cross-sections of (a) differential reflectivity and (b) nadir beam Doppler 
velocity for the Florida storm. S-band reflectivity contours are also shown. The 0°Clevel 
is indicated by dotted line. The thin dash line is the region with maximum attenuation. 
The thick dash line is the locus of the maximum updraft. 

Fig. 16. Vertical profiles of S-band reflectivity (solid), X-band nadir beam reflectivity 
(dashes), X-band nadir beam Doppler velocity (dash-dot), and differential reflectivity 
(dotted) for the Florida storm. The central plot is through the region of maximum path 
integrated attenuation as indicated by the dash line in Fig. 15, the left plot is 1 km left of 


it, and the right plot is 1 km to the right of it. 



Fig. 17. Similar to Fig. 15, for the Brazil storm. 


Fig. 18. Similar to Fig. 16, for the Brazil storm. 

Fig. 19. Normalized extinction cross-section of melting ice particles versus fraction of 
mass melted for selected melted diameters of 3, 4, 6, 7, 8, and 10 mm and three bulk 
densities of ice: a) 0.3 b) 0.5 and c) 0.92 g cm 3 . The ice particle is assumed to be 
spherical and melt is assumed to form a concentric sphere coat around it. The 
normalization is performed by the cross-section of the completely melted particle. 


Calculations are based on Mie theory. 
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Table 2. List of EDOP flight legs with P 
"*" are presented in the paper. 


Day 

Time 

Location 

* 980815 

2223-2237 

Florida 

980905 

2217-2221 

Florida 

990125 

2215-2223 

Brazil 

990125 

2242-2246 

Brazil 

* 990210 

1810-1815 

Brazil 

990212 

1955-2000 

Brazil 

990212 

2049-2055 

Brazil 


990217 


1846-1859 


Brazil 


larger than 20 dB. The legs indicated by 


Max. PIA (dB) 

cloud top (km) 

30 

18 

25 

16 

29 

16 

30 

16 

29 

15 

20 

13 

15 

14 


23 


13 








Table 4: Summary of inferred attenuation and other parameters for Brazil storm, 10 
February 1999. The values in the last four columns are for the height intervals in 3 rd 
column. 


Height (km) 
From top to 

Two-way 
PIA (dB) 

Height 

Interval 

(km) 

One-way 

attenuation 

(dB km ') 

Z s (dBZ) 

ZDR (dB) 

Z/k 

1.5 

21 
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23 
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1. 7-0.8 
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16 

2-3 

3.5 

-49 
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1.6-1. 7 

2.5X10 4 
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2.5X10 5 
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